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ABSTRACT 

Absorption between the rest-frame wavelengths of 973A and 1026A in quasar spectra arises from two 
sources (apart from occasional metals): one is due to Lya absorption by materials at a low redshift, and 
the other is from Ly/3 at a higher redshift. These two sources of absorption are to a good approximation 
uncorrelated because of their wide physical separation. Therefore, the two-point correlation of absorption 
in this region of quasar spectra neatly factorizes into two pieces: the Ly/3 correlation at high z, and the 
Lya correlation at low z. The latter can be independently measured from quasar spectra at lower 
redshifts using current techniques. A simple division then offers a way to statistically separate out the 
Ly/3 two-point correlation from the Lya correlation. Several applications of this technique are discussed. 
First, since the Ly/3 absorption cross-section is lower than Lya by about a factor of 5, the Ly/3 forest 
is a better probe of the intergalactic medium (IGM) at higher redshifts where Lya absorption is often 
saturated. Second, for the same reason, the Ly/3 forest allows a better measurement of the equation 
of state of the IGM at higher overdensities, yielding stronger constraints on its slope when used in 
conjunction with the Lya forest. Third, models of the Lya forest based on gravitational instability make 
unique predictions for the Ly/3 forest, which can be tested against observations. We briefly point out 
that feedback processes which affect higher density regions but leave low density structure intact may 
be better constrained by the Ly/3 forest. Finally, extending our technique to the higher Lyman series is 
in principle possible, but becomes increasingly difficult because of diminishing path lengths. 

Subject headings: cosmology: theory - intergalactic medium - large scale structure of universe; quasars 
- absorption lines 



1. INTRODUCTION 

A typical high redshift quasar spectrum is shown in Fig. 
1. Most studies of the intergalactic medium (IGM) focus 
on the part of the spectrum where the only kind of ab- 
sorption (aside from occasional absorption by metal sys- 
tems) is that by neutral hydrogen via Lyman-alpha (Lya) 
.i.e. the well-known Lya forest. This is the part of the 
spectrum that has a wavelength between (1 + z q )X^ and 

(1 + z g )A° , where z q is the quasar redshift. (Hereafter, 
we will use the symbols X° a = 1215.671, A° = 1025. 72A, 
and A^ = 972. 54A to denote the rest-frame Lya, Ly/3 
and Ly7 wavelengths.) Absorption blueward of (1 + z q )X'p 
is generally ignored because Lya absorption is entangled 
with absorption by the higher Lyman series (for excep- 
tions, see e.g. Press, Rybicki & Schneider 1993). For in- 
stance, absorption at a wavelength A between (1 + z g )A° 



IGM: 

u ~ 2c x (A° - Ag)/(A° + A° ) = 5.1 x 10 4 km/s (2) 

where c is the speed of light. The velocity separation u is 
much larger than the correlation scale at redshifts of a few: 
~ 100 km/s. This means that to good approximation, the 
Lya and Ly/3 optical depths that contribute to the total 
rtot. at a given observed wavelength are uncorrelated. This 
implies 



(g— Ttot.CA 1 ^— T t ot.(A 2 )^ 

= ( e -^(4)e-^(4))( e -^(4) e - r ' 3 (4)) , 



(3) 



and (1 + z q )X { l has two 



sources: 



-rtot. (A) 



(1) 



where r to t. is the net total optical depth observed at the 
wavelength A, and r a and t$ are the Lya and Ly/3 op- 
tical depths respectively, which arise from absorption by 
materials at two different redshifts: z a — [A/A^] — 1, and 
zp — [A/Ajg] — 1. We will refer to this part of the spectrum 
as the Lya + (3 region. 

The key observation of this paper is that z a and zp are 
associated with physically widely separated parts of the 



where (...) denotes ensemble averaging (or, operationally, 
averaging over lines of sight), and A 1 and A 2 refers to two 
different wavelengths that reside in the Lya + /3 region. 
The redshifts z^ and z 2 are the corresponding Lya red- 
shifts, and similarly, z^ and z^ are the corresponding Ly/3 
redshifts. 

The first equality in eq. (3), which states that the av- 
erage transmission in the Lya -I- (3 region factorizes into 
two parts, is implicitly assumed in the existing work that 
makes use of Ly/3 absorption (e.g. Cen & McDonald 2002, 
Fan et al. 2002, Lidz et al. 2002). 

The second equality goes one step further: it tells us 
that the two-point correlation in the same region also fac- 
torizes into two parts: the Lya correlation and the Ly/3 
correlation respectively. Since both the two-point correla- 
tion in e~ Ttat - and the two-point correlation in e~ Ta can be 
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Fig. 1. — Schematic diagram to show the relevant regions discussed in this paper. The region in which only Lya photons can be absorbed: 
((1 + z 9 )Ajj, (1 + 2 q )A^), where z q is the quasar redshift, is the famous Lya forest. The focus of this paper is on the Lyman a + (3 region, 

((1 + z q )\®, (1 + z q )X^). The spectrum is of the quasar Q2139-4434 at z q = 3.23, kindly provided by Arlin Crotts. Note that the rest-frame 
Lya, Ly/3 and hy-f wavelengths are: \° a = 1215.67A,A^ = 1025. 72A and A° = 972. 54A. The Lyman a emission line at A = 5144A is not fully 
shown. 



measured directly (the latter from separate lines of sight 
to quasars at lower redshifts) , eq. (3) tells us that we have 
a handle on the two-point correlation in e~ Tf > as well. 1 

The crucial point here is not to separate Lya and Ly/3 
absorption on an absorption-line by absorption-line basis, 
which is a challenging task. Rather, the strategy is to 
exploit the property of uncorrelated absorption to statis- 
tically separate out the two sources of absorption in the 
Lya + /3 region of quasar spectra. 

Why is the Ly/3 absorption interesting? Its utility lies 
in the smallness of the Ly/3 absorption cross-section: it is 
lower by a factor of 5.27 than the Lya cross-section. This 
means that for a given neutral hydrogen density, the as- 
sociated Ly/3 optical depth is a factor of 5.27 lower than 
the Lya optical depth. Therefore, Ly/3 absorption is more 
sensitive to structure at higher overdensities compared to 
Lya (because Lya goes saturated before Ly/3). As we ex- 
plain, this allows a better measurement of the equation of 
state of the IGM. Furthermore, by the same token, Ly/3 
absorption also offers a better hope for constraining the 
large scale structure at high redshifts (z > 5) where Lya 
absorption is often saturated. 

The organization of this paper is as follows. In §2, 
we demonstrate how the Ly/3 forest is more sensitive to 
high overdensities than the Lya forest, especially as far as 
the equation of state is concerned. In §3, we develop the 
method of statistically separating Ly/3 from Lya absorp- 
tion by rewriting eq. (3) in Fourier space. We demon- 
strate using a concrete example how the Ly/3 transmission 
power spectrum differentiates between different equations 
of state that have very similar Lya transmission power 
spectra. We end in §4 with a discussion of 1. how a better 
measurement of the equation of state might allow stronger 
constraints on the slope of the mass power spectrum, 2. 
how the Ly/3 transmission power spectrum provides a more 
precise test of the gravitational instability model of the 
forest, and yields more stringent constraints on feedback 
processes, and 3. corrections to eq. (3). 

1 Obviously, the same logic applies to the three-point correlation 
and so on, which wc will not discuss here. 



2. THE EQUATION OF STATE OF THE IGM: 
LYMAN-a VERSUS LYMAN-/3 

The photoionized IGM is well described by a temperature- 
density relation, or an effective equation of state, of the 
form: 

T = T A a (1) 

where T is the temperature, To is its value at mean density, 
a is the slope of the equation of state, and A = p/p, with 
p being the gas density, and p its mean. Such a relation 
appears to hold for A < 5 (Miralda-Escude & Rees 1994, 
Hui & Gnedin 1997). 

Current measurements from the Lya forest in the red- 
shift range of z <~ 2.4 — 4.0 yield much better constraints 
on To than on a (Ricotti, Gnedin & Shull 2000, Schayc 
et al. 1999, Bryan & Machacek 2000, McDonald et al. 
2001, Mciksin, Bryan & Machacek 2001, Zaldarriaga, Hui 
& Tegmark 2001 [ZHT01 hereafter]). The reason is quite 
simple to understand. Lya absorption is sensitive largely 
to A ~ 1 — 2, and the short lever arm does not al- 
low a precise measurement of the slope a. Present con- 
straints are consistent with the full physically plausible 
range a = — 0.6 (Hui & Gnedin 1997), according to 
ZHT01. 2 

To illustrate the range of overdensities to which the 
Lya forest is sensitive, we measure the transmission power 
spectrum from mock Lya spectra that are generated us- 
ing N-body simulations. The simulations used throughout 
this paper are 256 3 -grid, 128 3 -particle, P3M simulations 
of a Standard Cold-Dark-Matter model with a box size of 
16 Mpc/h. The simulations are dark-matter only simula- 
tions generated with the HYDRA code of Couchman et al. 
(1995). The baryon density is obtained by smoothing the 
dark matter density to mimic the effect of pressure forces 
(See ZHT01 for more details). These simulations have 

2 ZHT01 constrain a using measurements of the small scale flux 
power spectrum. McDonald et al. (2001), using a line-fitting 
method, find a tighter constraint on a at z = 2.4, requiring a > 0.38 
at a lcr-confidence level. Schaye et al. (2000) obtain still tighter 
constraints also using a line— fitting method. The difference between 
the Schaye ct al. (1999) error bars and the McDonald et al. (2001) 
error bars is due to different line selection criteria (See McDonald ct 
al. (2001) for a discussion of this). 
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Fig. 2. — Illustration of how Lya and Ly/3 absorption probes different overdensities. Left: Differential Lyo transmission power spectra for 
a set of artificial temperature-density relations defined by eq. (1) for A c ranging from to 3.6 in steps of 0.4 (bottom to top, as seen around 
k = 0.1 s/km). What is shown in each case is not the absolute power spectrum, but rather the logarithmic difference from a baseline power 
spectrum denned by that of A c = (this is why the line corresponding to A c = is exactly zero). Right: Differential Ly/3 transmission power 
spectra for A c ranging from to 6.0 in steps of 0.4. The redshift in both panels is z = 3.0. [See the electronic version of the journal for a 
color version of this figure] 



lower resolution than recommended by McDonald (2003) 
and Meiksin & White (2003) to achieve convergence in 
measuring the flux power spectrum. We expect, however, 
that our present simulations are adequate to qualitatively 
illustrate the benefits of measuring the Ly/3 flux power 
spectrum. 

In Figure 2 {left panel) we show the Lya transmission 
power spectrum at z = 3.0 for a set of artificial equations 
of state parameterized in the form (following ZHT01): 

T _ / T for A < A c , . 

1 ~ \ 2T for A > A c [Z> 

where A c defines the density threshold where the gas 
temperature T has a jump and To = 1.2 x 10 4 K. The left 
panel shows the Lya transmission power spectrum for A c 
varying from 0.0 to 3.2 in steps of 0.4. We have subtracted 
from each power spectrum a baseline power spectrum cor- 
responding to that of A c = 0.0. In comparing power spec- 
tra with different A c , we fix all of the other IGM modeling 
parameters. 

The Lya transmission (or flux) power spectrum Pt a (k) 
is defined as: 

€f,a(u) = (6f, a (u 1 )6 f , a (u 1 + u)) , (3) 

where £,f, a (u) is the two-point correlation function at a 
velocity separation of u, the transmission power spec- 
trum P/>(fc) is its Fourier counterpart as a function of 



wavenumber k, and Sf_ a is the fluctuation of the Lya trans- 
mission about its mean: 



Note that the definition of £/ jC , differs from the kind of 
two-point correlation in eq. (3) in that £y iQ involves the 
correlation of which has zero mean, and is properly 
normalized. We will sometimes refer to £/ as the normal- 
ized two-point function. 

Figure. 2 (left panel) shows that the Lya transmission 
power spectrum at high k (~ 0.1 s/km) increases, as A c 
increases from (or, in other words, according to eq. (1), 
as the temperature drops). This is consistent with the 
picture that thermal broadening smooths the flux field, 
and suppresses small scale power, a fact that is used to 
constrain the thermal state of the IGM from observations 
(ZHT01). The interesting point of Figure 2 (left panel) is 
that changes in the Lya transmission power spectrum are 
most pronounced in a narrow range of A c 's, from about 1 
to 2. For smaller or larger A c 's, the Lya power spectrum 
does not change in a substantial way. This is fundamen- 
tally why current observations from the Lya forest yield 
useful constraints on the temperature around mean den- 
sity (or 1 — 2 times the mean density), but not on the 
slope of the equation of state a (eq. (1) ). A stronger con- 
straint can be obtained if one has the means to probe the 
temperature of the IGM over a larger range of densities. 

The Ly/3 transmission power spectrum provides the req- 
uisite probe of the IGM over a large range of densities, as 
illustrated by Figure 2 (right panel). This panel is similar 
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to the left panel, except that it shows the Ly/3 transmis- 
sion power spectrum rather than the Lya power spectrum. 
The Ly/3 transmission power is defined by (similar to eq. 
(3), and (4) ): 



OA u ) = ( s fA uX ) 5 fA ul + u )) > 



(5) 



(6) 



In the right panel, we let A c vary from 0.0 at the bottom 
to 6.0 at the top. Clearly, the Ly/3 transmission power 
spectrum is sensitive to higher densities: it continues to 
vary from A c of 1 to about 4. 3 

Note that the experiment shown in Fig. 2. is done at 
z = 3. If one moves to higher redshifts, the A one is 
sensitive to shift to lower values for both Lya and Ly/3. 
This is simply because the mean density of the universe 
is higher, and lower overdensities give rise to the same 
amount of absorption as at z — 3. For instance, we have 
repeated the above experiment at z — 3.74, and find that 
the range of sensitive A shifts downward by about 0.5 for 
both Lya and Ly/3. 

Fig. 2. is only meant to illustrate the range of densities 
to which Lya and Ly/3 absorption are sensitive. Does the 
Ly/3 transmission power spectrum differentiate between re- 
alistic equations of state that are otherwise difficult to dis- 
entangle using the Lya power spectrum alone? This is ad- 
dressed in Fig. 3. At the top of the figure, we show three 
different Lya transmission power spectra, each assuming 
a different equation of state slope (a in eq. (1) ). The 
other IGM model parameters (see ZHTOI) are adjusted 
slightly (well within observational uncertainties) 4 in each 
case to give very similar P/ jQ (fc)'s. The Ly/3 transmis- 
sion power spectra for these three cases show more differ- 
ences at high k, suggesting that observational constraints 
on Pfj(k) might be useful. Our next task is to describe 
how to tease out Pf t p(k) from the Lya + /3 region of quasar 
spectra. 

3. CONSTRAINING THE LY/3 POWER SPECTRUM 

As discussed in §1, the coincident Lya and Ly/3 absorp- 
tion in the Lya + (3 region of a quasar spectrum are to a 
good approximation uncorrelated. This gives rise to a sim- 
ple factorization of the (un-normalized) two-point function 
(eq. (3) ). The normalized two-point function obeys 

O,tot» = €fA u ) + €fA u ) + tfA u )€fA' u ) C 1 ) 

where £/,tot. as a function of velocity separation u is de- 
fined in a similarly to eq. (3) and (5) (i.e. ^/,tot.(w) = 
(*/,tot.(ui)*/,tot.(«i+«)),with 6/,tot. = [e~T**-/(e-T**-)]- 

3 At sufficiently high overdensities, the equation of state given in 
eq. (1) would no longer be a good description of the temperature- 
density relation, because large— scale shock heating becomes impor- 
tant, which introduces a significant scatter to the temperature. The 
temperature-density relation for overdensities of A up to ~ 5 can 
probably still be described by eq. (1) (Hui & Gnedin 1997). 

4 In the notation of ZHT01, model A/B has (k F ,a,T ,f) = 
(36,0.0,310,0.512) / (39,0.4,305,0.507). Here k F is in units of h 
Mpc -1 and Tq is in units of (km/s) 2 . 
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Fig. 3. — Transmission power spectra for Lya (top three curves) 
and Ly/3 (bottom three curves) for three models that have different 
equations of state. The models are chosen so that the Lya power 
spectra are very similar. This is at z = 3.74. [See the electronic 
version of the Journal for a color version of this figure] 



1, and Ttot. being the total Lya + Ly/3 optical depth). The 
velocity separation u is related to the quantities A 1 , A 2 , z\ 
z 2 , Zp and z"^ in eq. (3) by: 



A 2 - A 1 zl 
c = = c- 



A 



l + z n 



= c ^i (2 ) 



1 



Z0 



where A is the mean observed wavelength in the Lya + (3 
region, and A = (1 + z Q )A" = (1 + Z/s)Ajg. 
The Fourier counterpart of eq. (1) is 

P f ,,ot.(k) = P f , a (k) + P f Ak) (3) 

+ / ^PfAk-k')P f Ak') 

where P/,tot. is the power spectrum of the total transmis- 
sion in the Lya + /3 region. We will often refer to f/,tot. as 
the total power. Note that implicit in the above expres- 
sion is that Pf >a and Pf,p are at different mean redshifts: 
z a and zp. 

One can directly measure both P/,tot. and Pf, a from 
observational data. To constrain the Ly/3 power spectrum 
-P/,/3, it is important that quasars at different redshifts 
are employed to measure -P/,tot. and Pf, a - the former 
from higher redshift quasars, and the latter from lower 
redshift quasars whose Lya-only region overlaps in wave- 
lengths with the Lya + /3 region of the former. 

In principle, once Pf, tot. and Pf >a are given, eq. (3) can 
be inverted to obtain -P/,/3 if one thinks of it as a linear 
vector equation: 



Pf.tot. - Pf,a = M • Pf i/3 



(4) 
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Fig. 4. — Fractional differences (squared) between the two models 
A and B (with the slope of the equation of state a = 0.0 and a = 
0.4 respectively; see Fig. 2.) in the total and Lya power spectra. 
[See the electronic version of the Journal for a color version of this 
figure] 

where M is a matrix whose components are 

dk ■ 

M(k u kj) = S i:j + -lp f fa - kj) . (5) 
Ztt 

In other words, from eq. (4), one can in principle ob- 
tain: 

P fl/3 - M- 1 • (P f , tot . - P f>a ) (6) 

While such an inversion is useful for visually inspecting the 
Ly/3 power spectrum, in practice it can be noisy and one 
is likely better off focusing on the total power but keeping 
in mind that the Lya contribution is known. 

We therefore will not pursue the path of inversion here. 
Instead, we will be content with posing the question: how 
different is the total observable power P/,tot. for models 
that are quite degenerate in their Lya power, for instance 
those shown in Fig. 3? 

The answer is provided by Fig. 4. Shown here is 
the fractional difference (squared) in the total power (top 
curve) between the two models labeled A and B (a = 
and a = 0.4 respectively) in Fig. 3. The total power here 
is computed using eq. (3) by taking Pf^ at zp — 3.74, and 
at — 3.0. 

For comparison, we show in the same figure the frac- 
tional difference in the Lya power (bottom curve) between 
the same two models at the same redshift as the Ly/3 
power, i.e. z = 3.74. Comparing the two curves tells us 
how much more the total power (due to the information 
from Ly/3 absorption) can constrain the equation of state 
at z = 3.74 compared to the Lya power spectrum alone. 



We assume implicitly here that Pt a at z a = 3.0 is well 
known from data at low redshifts. 

From Fig. 4, we can infer the signal-to-noise ratio (S/N) 
with which one can distinguish between the two different 
equations of state (models A and B). The ( S/N) 2 is given 

by 

(S/N) 2 = J2 i P f( k ) ~ P ?( k )} 2 / M*0 2 (7) 

k 

where P^ and P? are the power spectra for model A 
and B respectively (Pf here can stand for either Pf. a or 
-P/,tot.)> and ap(k) 2 is the variance in power (here taken 
to be that for model B). The summation extends to all 
modes with positive k values, up to fc max = 0.12 s/km. 
The power at higher k is expected to be contaminated sig- 
nificantly by metal lines (McDonald et al. 2000) and so 
the sum is truncated at fc max . The variance o~p(k) 2 equals 
(P?(k) + shot) 2 , if one approximates the fluctuations as 
Gaussian random (see Hui et al. 2001). It is likely inac- 
curate to assume Gaussian errors for some of the high k 
modes considered here. A more accurate error estimate is, 
however, unlikely to change our main point that the Ly/3 
flux power spectrum distinguishes between models with 
different temperature-density relations, a, that have iden- 
tical Lya flux power spectra. 

The shot-noise term 'shot' is approximately equal to 
(Au/(e~ T ))(s/n)~ 2 , where Au is the size of a pixel in 
km/s, and s/n is the signal-to-noise ratio per pixel of the 
data (different from S/N defined above). High quality 
spectra (e.g. Kim et al. 1997, Kirkman & Tytler 1997, 
Simcoe, Sargent & Rauch 2002, Kim et al. 2002) often 
have a sufficiently low shot-noise that shot-noise can be 
ignored altogether, in which case, (S/N) 2 — Ylk(Pf(k) ~ 
Pf (k)) 2 / Pf (k) 2 , which is the area under the curve shown 
in Fig. (4). 

We find that using the Lya power spectrum at a mean 
redshift of z — 3.74 gives a discriminating power of 
(S/N) 2 = 0.8, while using the total power from the Lya+/3 
region (corresponding to Ly/3 at the same z = 3.74) gives 
(S/N) 2 = 7.6. This is the (S/N) 2 for one line of sight, as- 
suming that the Lya and Lya + /3 regions have a length of 
5.1xl0 4 and5.1xl0 4 (A°-A°)/(A^-A°) = 1.4xl0 4 km/s 
respectively. The discriminating power (S/N) 2 scales lin- 
early with the number of line of sights used. Note that 
this estimate assumes that one has enough resolution to 
measure the power spectra at high k (up to k = 0.12 
s/km) i.e. a resolution of Full-Width-at-Half-Maximum 
(FWHM) = 8 km/s or better. In summary, using the Ly/3 
forest boosts the discriminating power (S/N) between an 
equation of state of a = and an equation of state a = 0.4 
by a factor of ~ 3. If only modes up to k = 0.1 s/km are 
included then this number goes down to 2.5. These num- 
bers are typical for models very close in Lya. 

To completely quantify how much better one can con- 
strain the equation of state, a, using the Ly/3 flux power 
spectrum in addition to the Lya flux power spectrum, a 
more rigorous approach is necessary. One should run a 
large grid of simulated models, simulating both P/,tot. and 

5 In order to calculate the total power, we use the same P/ lQ at 
z a = 3.0 in model A and model B - the difference in the total power 
of these models comes only from differences in their Ly/3 power. 
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Pf, a , compare with mock data, and marginalize over all of 
the other modeling parameters to obtain the reduced like- 
lihood function for a. A comparison between the result- 
ing likelihood function formed using Pf _ a alone, with that 
formed using P/.tot. plus P/, Q) quantifies how much one 
gains using the Ly/3 forest. Here we are content to illus- 
trate that the Ly/3 flux power spectrum adds extra infor- 
mation on a that is not available from the Lya flux power 
spectrum alone (Figure 3), without completely quantifying 
how much tighter the resulting constraints on a are. 

4. DISCUSSION 

The discussion §3 suggests that the Ly/3 forest can in- 
deed be beneficially used in conjunction with the Lya for- 
est to help discriminate between different models, espe- 
cially between different equations of state. The reason 
for its utility lies in its sensitivity to higher overdensities. 
Better measurements of the equation of state are useful in 
constraining the rcionization history of the universe (The- 
uns et al. 2002, Hui & Haiman 2003). Since spectra of 
quasars at sufficiently high redshifts often extend well into 
the Ly/3 forest, there is no reason not to exploit this part of 
the spectrum to increase the scientific return. An alterna- 
tive possibility for constraining the equation of state is to 
use the Lya flux bispectrum in conjunction with the flux 
power spectrum. Referring to Table 3 of Mandelbaum et 
al. (2003), it appears that combining the flux bispectrum 
and power spectrum yields tighter constraints on the tem- 
perature density relation than using the power spectrum 
alone. 

Zaldarriaga, Scoccimarro & Hui (2001) pointed out that 
from the Lya power spectrum alone, there is a near- 
degeneracy between models that trade off variations in the 
mass power spectral index n s with variations in the equa- 
tion of state. Our finding that the Ly/3 forest can place 
stronger constraints on the equation of state suggests that 
one might be able to break this degeneracy using the Ly/3 
power spectrum. 

We have carried out an experiment similar to that in Fig. 
4, except we replace models A and B with the following 
two models: one has n s = 0.7, a = 0.0 and Tq = 2.1 x 10 4 
K, and the other has n s = 1.1, a = 0.6 and To = 1.6 x 10 4 
K (see eq. (1) for definitions of a and T ); these two 
models have quite similar Lya power spectra. The total 
power distinguishes between the models at a level similar 
to that shown in Fig. 4. However a perhaps more readily 
realizable option is to break the degeneracy by measuring 
the Lya power spectrum as accurately as possible on large 
scales, k ~ 0.001 — 0.01 s/km. In this range the models 
are different in the Lya power spectrum, albeit at a level 
that is too small to distinguish with existing data. Quasar 
spectra from the Sloan Digital Sky Survey are well suited 
for this. 

An interesting use of the Ly/3 forest is to search for signs 
of feedback processes in the IGM. Recent interest in the 
Lya forest as a cosmological probe relies on a framework in 
which fluctuations in the forest arise naturally from grav- 
itational instability (e.g. Bi, Borner, & Chu 1992, Cen 
et al. 1994, Zhang et al. 1995, Reisenegger & Miralda- 
Escude 1995, Hernquist et al. 1996, Miralda-Escude et al. 

1996, Muecket et al. 1996, Bi & Davidsen 1997, Bond & 
Wadsley 1997, Rauch et al. 1997, Hui, Gnedin & Zhang 

1997, Croft et al. 1998, Theuns et al. 1999, Nusser & 



Haehnclt 2000, McDonald et al. 2000, White & Croft 
2000, Meiksin et al. 2001, Pichon et al. 2001, Croft et 
al. 2002, Gnedin & Hamilton 2002, Viel et al. 2002). An 
important assumption behind this picture is that feedback 
processes, such as galactic winds (Adelberger et al. 2003), 
do not significantly disturb the IGM. The good agreement 
between observations (particularly the Lya transmission 
power spectrum) and the gravitational instability model is 
often used as an argument that feedback processes, while 
inevitably present, do not affect large volumes of the IGM. 
A reasonable expectation is that they preferentially af- 
fect higher density regions. If so, the Ly/3 forest offers 
a better hope of testing for the presence of such feed- 
back processes. Particularly interesting is the fact that 
once the gravitational instability model parameters (such 
as the mass power spectrum, cosmology, etc) are tuned to 
match observations of the Lya forest, there are definite 
predictions for the correlations observed in the Ly/3 for- 
est, since Lya and Ly/3 optical depths are simply related 
by a rescaling in the cross-section. As we have discussed, 
changing the slope of the equation of state (a) does seem 
to modify the Ly/3 power spectrum while leaving the Lya 
power spectrum relatively unchanged. Feedback processes 
might behave in the same way i.e. affecting the Ly/3 forest 
more than the Lya forest, except that it is unlikely that 
their effects can be mimicked by simply varying a. Galac- 
tic winds for instance change the density structure of the 
IGM, by creating evacuated shells around galaxies. They 
might also introduce a larger than expected scatter in the 
IGM temperature at high overdensities. 

It is important to reiterate the method we advocate is 
not to decipher the Ly/3 forest on an absorption line-by- 
absorption-line basis. Rather, the strategy is to statisti- 
cally detect the presence of Ly/3 correlation, exploiting the 
fact that the Lya and Ly/3 absorptions that fall within the 
Lya + (3 region of a quasar spectrum are uncorrelated (eq. 
(3) ). A natural question is: how good an approxima- 
tion is it? The fractional correction to the first equality 
of eq. (3) is the two-point correlation between Lya and 
Ly/3 transmission fluctuations at a velocity separation of 
u <~ 5.1 x 10 4 km/s (eq. (2) ). It is safe to assume that 
this correlation is weaker than the Lya two-point corre- 
lation (£/.a; see eq. (3) ), since we know from experi- 
ence that Ly/3 absorption weakens the correlation (com- 
pare the Lya and Ly/3 curves in Fig. 3). The observed 
two-point correlation function in the Lya forest has only 
been reliably measured at velocity separations of u <; 1700 
km/s (McDonald et al. 2000). In order to estimate the 
two-point correlation function at velocity separations of 
ii ~ 5.1 x 10 4 km/s, we extrapolate from our simulation 
measurements assuming linear biasing (Scherrer & Wein- 
berg 1998, McDonald et al. 2000, Lidz et al. 2002). The 
resulting estimate is £/>(" = 5.1 x 10 4 km/s) <; 10~ 5 . 

Therefore, any correction to the first part of equation (3), 
(g-Ttot.wj _ ( e -T a (z a )^ e -/3 a (zf,)} ; must b e very sma n : 

<; 10~ 5 . Such a correction is even smaller at lower red- 
shifts where the two-point correlation is weaker. We can 
similarly estimate the corrections to the factorization of 
the two-point function (second part of eq. (3), or equiva- 
lcntly, eq. (1) ). The fractional error we make in the two 
point function is roughly ~ 1£,f,ap{u a p,z a p)/£,f,p{u,zp). 
Here £/, a /3 refers to the two point correlation between the 
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Ly/3 absorber that absorbs at observed wavelength A 2 and 
the Lya absorber that absorbs at A 1 , z a p is the mean red- 
shift between these absorbers, and u a p is their velocity 
separation. The fractional error depends on the relative 
size of u and u a p. Here we estimate the error when the 
Ly/3 absorbers are separated by one correlation length, 
u <~ 100 km/s, and situated at zp = 3.74. In this case 
the Lya absorber is separated from the Ly/3 absorber by 
u a fj ~ 5.1 x 10 4 km/s. A conservative error estimate 
then comes from taking £f,p(u,zp) ~ £f ta (u,zp)/5, (see 
Fig. 3), and £f, a p(u a p,z a p) ~ £f, a (u a p,zp)/\/E. From 
McDonald et al. (2000), the observed two-point func- 
tion at a separation of one correlation length is £/, Q (it ~ 
100 km/s, z ~ 4) ~ 0.2, which shows that the fractional 
error is £ (10/v / 5)(10 _5 /0.2) ~ 0.03%. For very widely 
spaced pixels in the quasar spectrum, the fractional error 
is larger, but these pixels have a negligible correlation. 

Two issues are worth further exploration. First, one 
must take care in masking out the intervening metal ab- 
sorption lines in the Lya + j3 region of the spectrum, just 
as one normally does in the Lya forest. Metal lines that 
cannot be easily masked out, such as OVI, can be suitably 
divided out using a method similar to the one mentioned 
in this paper i.e. use lines of sight where the red-side of 
Lya coincides in wavelengths with the Ly/3 forest of in- 
terest; the two-point correlation of metal absorption from 
these lines of sight can be used to take out the metal con- 
tamination, much as we remove the Lya contamination to 
the Ly/3 forest (a similar technique for 'cleaning' the Lya 
forest has been developed by P. McDonald and U. Seljak 
[2003, private communication]). Second, it is interesting to 
explore whether the ideas presented here can be extended 
to the higher Lyman series. While this is in principle pos- 
sible, it is likely that the combination of diminishing path 
lengths, and the increasing entanglement of different kinds 
of absorption (i.e. the coexistence of Lya, /3, 7, and so on), 
makes it difficult to exploit the factorization of correlations 
in practice. 
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